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Abstract 

We review the construction of exactly solvable lattice models whose 
continuum limits are N = 2 supersymmetric models. Both critical and 
off-critical models are discussed. The approach we take is to first find 
lattice models with natural topological sectors, and then identify the 
continuum limits of these sectors with topologically twisted N = 2 
supersymmetric field theories. Prom this, we then describe how to 
recover the complete lattice versions of the N = 2 supersymmetric 
field theories. We discuss a number of simple physical examples and 
we describe how to construct a broad class of models. We also give 
a brief review of the scattering matrices for the excitations of these 
models. 



USC-93/026 
hep-th/931lT38 



November, 1993 



1 Introduction 



N = 2 supersymmetric theories in two dimensions are in some ways much simpler than 
their non-supersymmetric kin. This is essentially because the N = 2 supersymmetry 
imples the presence of a topological sector for which semi-classical analysis yields exact 
quantum results. (In the language of supersymmetry, the "F-terms" are not renormal- 
ized.) This does not mean that the model is semi-classically rigid, but only that a key 
(topological) subsector of the theory is semi-classically determined: the complete theory 
is very rich and has all of the complexity of a non-supersymmetric theory. The topo- 
logical subsector has thus provided a "bridgehead" from which many of the non-trivial 
quantum aspects of the model can be explored, and usually with greater facility, and 
often in more detail than is possible for non-supersymmetric theories. (For reviews, see 
[|^, 0, 0.) For example, N = 2 Landau- Ginzburg theories have a superpotential that is 
exact. For N = 2 superconformal models this means that the operator algebra and the 
anomalous dimensions (conformal weights) of the Landau-Ginzburg fields are trivially 
computable 0, ^. On the other hand, it is still somewhat unclear as to how much of 
the rest of the theory, and in particular, how much of the complete operator content is 
determined by the Landau-Ginzburg potential. Indicative results are known: For ex- 
ample, the ADE classification of modular invariants of = 2 minimal models collapses 
to the ADE classification of modality zero singularities 0, ^ . Ramond sector characters 
can be obtained from the Landau-Ginzburg potential by computing the elliptic genus 

i I, i m- 

For N = 2 supersymmetric quantum integrable field theories (QIFT'sQ) the Landau- 



Ginzburg description leads to transparent analysis of the soliton structure [|TI], |T^, [T^. 
In A^ = 2 QIFT's, and even non-integrable N = 2 QFT's, differential equations can also 
be obtained for some of the scaling functions [Q. For the N = 2 QIFT's these scaling 
functions can also be obtained from the thermodynamic Bethe ansatz, but instead of 
differential equations, one obtains complicated integral equations |T3 



Over the last decade it has also become evident that many of the structures of 
quantum integrable field theories have analogues in exactly solvable lattice models (see. 



for example ||16| , \i7l [181 , p!9| , [20||). Using this technology, lattice models have been 
constructed in which the continuum limits give rise to many of the non-supersymmetric 
conformal field theories. It is therefore natural to expect that there should be exactly 
solvable lattice models whose continuum limits are N = 2 supersymmetric QIFT's. We 
shall henceforth refer to such lattice models as A^ = 2 lattice models^. It is also to 

^In this review, conformal field theory will be abbreviated as CFT, quantum field theory 
as QFT, and quantum integrable field theory as QIFT. 

^This nomenclature is somewhat misleading in that it suggests that the supersymmetry 
is realized on the lattice, whereas it still remains unclear whether this can be done in most 
of the N — 2 lattice models thus far constructed. 
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be hoped that such lattice regularizations may have especially simple features. If so, 
they may in turn help us to understand some mysterious issues like the appearance of 
branching functions in local height probabilities ||21|. There are also practical reasons 
for constructing N = 2 lattice models. The simplicity of = 2 QFT's has led to 
more progress than in any other kind of field theory. It is especially attractive to try to 
compare some of the advanced results with real or computer experiments. To do so, the 
construction of a "physically reasonable" lattice model whose continuum limit is the 
N = 2 QFT of interest is a natural way to proceed. We also hope that the N = 2 lattice 
models will be "more universal" than their non-supersymmetric counterparts. Indeed, 
consider a general non-supersymmetric model, and suppose that we do not impose 
the requirement of exact solvability. Since there are generically very many relevant 
operators in such a model, the corresponding coupling constants would need to be very 
finely tuned in order to find a particular second order phase transition and hence a 
particular conformal field theory. Thus one of the side-effects of exact solvability is to 
automatically make such a fine-tuning of the couplings. The N = 2 lattice models have 
a further special, and robust, identifying feature: the presence of a topological sector. 
Thus, not only can one construct them by imposing exact solvability, but one can also 
identify them by virtue of their topological sector. In fact, we will argue in this review 
that any "trivial" statistical mechanics model will probably give rise to a non-trivial 
N = 2 lattice model, and thus the N = 2 models will be of relevance in a large variety of 
situations. For instance the polymer and percolation problems, the statistics of Bloch 
walls in high temperature Ising model, Brownian motion, are all described by = 2 
supersymmetric theories. 

The investigation of lattice models associated with supersymmetric QFT's has a 
rather long history. To our knowledge^ it started with the study of = 1 supersym- 



metry in the tricritical Ising model [^, ^ ^ and the fully frustrated XY-model 
(the latter related to Josephson junctions arrays). 

Then = 1 supersymmetry together with N = 2 supersymmetry were observed 
at some special points of the Ashkin- Teller model and the six-vertex (or XXZ) model 
[p^ , |27| , ^] . By analysis of torus partition functions it was easy to identify more gener- 
ally special points of the higher spin Uq{SU{2)) vertex models p9|, [S^] that were N = 2 
supersymmetric. A physical explanation of this came later [^. In the simplest case of 
the six-vertex model it leads to A^ = 2 supersymmetry in the percolation problem. Sim- 
ilarly, by analysis of the related Izergin-Korepin (0(n)) model |3^, |33[, another family of 
lattice models with N = 2 supersymmetry was found (or conjectured). The simplest of 
this family leads to A^ = 2 supersymmetry in the self- avoiding walk (polymer) problem. 
■^We apologize for any reference missing in this short history. 
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This turns out to be a very favorable case for comparison with real and numerical ex- 



periments 1^ . More complete study of the SU (2) based models was carried out in |35 



Recently, it has been shown how one can go beyond models based upon SU (2). That is, 
N = 2 lattice models based on any simply-laced Lie algebra have been constructed by 
using "partial restriction" of modified solid-on-solid (SOS) models, or equivalently by 
twisting and performing partial quantum group truncation of the corresponding vertex 
models. The continuum limits of these lattice models are the N = 2 superconformal 
coset models based on hermitian symmetric space |^ |3^ . 

At the present time, there is a respectable number of known N = 2 lattice models. 
They resemble in some respects the corresponding QFT, although the structures are 
not yet as closely linked as one would wish. They have already led to comparison with 
real and numerical experiments, and very good agreement has been found. The purpose 
of this review is to describe in some detail the current state of knowledge and to point 
to new directions of research. 

The first section is rather qualitative and collects observations about the expected 
structure of iV = 2 lattice models. A variety of simple examples is worked out, starting 
from simple geometrical ideas and introducing the fundamental tools that will be used 
extensively later. One of the basic themes underlying the N = 2 lattice constructions is 
the role of the "trivial" topological sector. This will be used directly in later sections to 
construct and analyse the N = 2 lattice models based on general Lie algebras. In the 
fourth section we will briefly discuss the scattering matrices of excitations in the N = 2 
lattice models and their continuum limits. In section five we discuss further physically 
interesting aspects of = 2 lattice models, and in the last section we conclude by 
summarizing what we believe are the important open problems in the subject. 



2 Features of = 2 lattice models. 

Our purpose in this section is to study a number of closely related and simple models 
that exhibit the basic features of the N = 2 lattice models. Our aim is to try to give 
some intuitive understanding of A^ = 2 lattice models, to show how one can analyse 
their various features, and ultimately to evolve a strategy for finding such models. 

2.1 The clues from the topological sector 

If the theory has N = 2 supersymmetry then, as mentionned earlier, there is a topo- 
logical sector, with the following characteristics. The topological sector of the N = 2 
QFT consists of operators that are annihilated by two of the four supercharges. The 
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topological physical^ states consist of only the ground states of the non-topological 
N = 2 lattice model. The operators of the topological model are order parameters of 
the non-topological sector, and the correlation functions of these operators are constant 
in the topological sector. 

The presence of some sort of "topological" sector is one of the main attributes 
of an = 2 lattice model, and because of the links between the topological and 
non-topological sectors in the continuum, the identification of the lattice topological 
sector proves an important tool in the construction of the complete N = 2 lattice 
model. More precisely, the idea is to study representations of the lattice algebras and 
lattice symmetries (such as Temperly-Lieb and quantum group representaions). For 
certain choices of parameters some of these will be very simple or trivial, and the 
corresponding lattice models will be "topological." The idea is to define the topological 
model precisely, fix the parameters to their topological values, but then modify the 
choice of representations or modify the choice of lattice variables in such a manner that 
the N = 2 lattice model emerges from its topological limit. We will describe two "dual" 
methods for accomplishing this. The first relies on the fact that a lattice model can have 
very different properties when it is described by different sets of observables that are 
not mutually local. One thus obtains some of the N = 2 lattice models by considering a 
"trivial" model where the given lattice variables do not interact, and then one makes a 
"change of variables" to new set of geometrical observables that have highly non-trivial 
properties. The non-interacting observables describe the topological model, while the 
new observables describe the non-trivial N = 2 lattice model. The second method is 
to consider the restriction process in solid-on-solid (SOS) models (or quantum group 
truncation in the equivalent vertex formulation), and start with a restriction that results 
in a "frozen" or completely rigid topological model. One then "gets something from 
nothing" by relaxing the restrictions imposed on the SOS model, while at the same time 
one "topologically untwists" the Boltzmann weights to avoid singularities in the transfer 
matrix. It turns out that the first method is a priori a little more intuitive, while the 
second method is much easier to generalize. We therefore start by describing a simple 
example of the first procedure and then trace a connection to the second method. 

2.2 Examples of lattice topological sectors 
2.2.1 Percolation: an example of a geometrical model 

Consider a square lattice and put on every edge a variable a that can take two possible 
values, or 1. Suppose that the edge variables do not interact with each other, but put 

''^The definition of "physical" for the lattice model can be quite different that of the 
field theory. 
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the whole system in a magnetic field. Since the spins do not interact, all correlation 
functions are trivially constant (or, at least are delta functions). This appears to be 
a completely trivial model and the partition function is Z = il + e j , where Ne 
is the number of edges. We now pass to a new set of observables (called geometrical 
observables in what follows): Call edges with variable a — "empty" and edges with 
variable a — 1 "occupied" (we represent occupied edges with heavy lines in figure 1) 
and consider the geometrical properties of clusters, that is, of the connected sets of 
occupied edges. These cluster observables are non-local with respect to the original a 
variables, and the "trivial" model now describes the bond percolation problem. Since 
the sum over variables cr = 0, 1 is unconstrained, there is a probability of p — ^^JIh 
of having an occupied edge, and a probability of p = ^_^_1-h of having an empty edge. 
The correlation functions of new variables, such as the probability that two vertices 
belong to the same cluster, are non-trivial. The original trivial model can, of course, 
be recovered by deciding to study only local questions like "is this edge occupied" 
and forgetting about non local questions like "are these two edges part of the same 
cluster". As we will see, the distinction between these two questions is clearer from a 
representation theory point of view, and easy to implement, for instance, in a transfer 
matrix formalism. 

To determine the supersymmetry point of this model we consider it upon a torus. 
The new observables can be given various boundary conditions: for instance, one can 
sum over non-contractible clusters giving each some weight, or fugacity. The presence of 
unbroken supersymmetry means that, in the thermodynamic limit, the free energy per 
edge cannot depend upon the boundary conditions on the torus. Intuitively, if p > 1/2 
or p < 1/2, one will find large clusters of occupied or unoccupied edges. Thus non- 
contractible clusters will make significant contributions to the free energy, and thus the 
latter will depend upon the boundary conditions. However, when one has p = 1/2, or 
H — 0, one is at the percolation point, a critical point where the geometrical correlation 
functions decay algebraically and there is a vanishing fraction of big clusters of occupied 
or unoccupied edges. The partition function is Z = 2^^. This lattice partition function 
is, as usual, defined up to a non-universal term of the form e-^^^. In the continuum, the 
partition function simply counts the number of "different observables" of the model. It 
is rather reasonable that there is only one such observable (beside the identity): the 
variable a. One thus expects Z = 2. This is also consistent with the identification 
of this percolation problem as an = 2 lattice model. As we will see, the foregoing 
lattice partition function can be interpreted as the Witten index of the N = 2 lattice 
model. We will soon present rather more direct evidence that for p = 1/2 this is indeed 
an A^ = 2 lattice model. 
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Figure 1: a) A configuration of "up" (black dot) or "down" (wliite dot) spins 
on the edges of the square lattice, b) The corresponding configuration of 
"occupied" or "empty" edges. 



2.2.2 The dual point of view: the "frozen" Potts model 

The foregoing non-trivial bond percolation problem can also be described using the 
Q-state Potts model in the limit Q — > 1. The general Q-state Potts model is defined 
by putting a spin variable, a, on each vertex of the square lattice. These spin variables 
take values in 1, . . . ,Q, and the lattice partition function is 



^Potts — TT 

{cr} <i;j> 



S{ai,aj) 



(1) 



where the sum is over all spins and the product is over is over nearest neighbours. 
The connection with the percolation problem becomes evident upon making a high 
temperature expansion of the Potts model: The partition function can be written 



^Potts 



tr 



( n {l + (e^-l)<^(<T„a,)} 

\<i;j> , 



(2) 
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which can then be expanded graphically by connecting all vertices with the same spin 
a. One then finds that 

Z = ^ (e^ - 1)^^ , (3) 

graphs 

where Nb is the number of bonds in the graph, and Nc is the number of clusters. This 
high temperature expansion coincides with the bond percolation problem only if the 
Q'^'^ term is equal to one, i.e. if and only ii Q = 1. One must also make the following 
identification: 

- 1 = e-^ . (4) 

However, if one considers the Potts model with Q = 1, one finds that it is a trivial, 
completely frozen, "topological" model. The partition function collapses to ^potts = 
qKNe_ ^]^g bond percolation problem should really be viewed as the Q ^ 1 limit of 
the Potts model. For example, the derivatives of ^potts at Q = 1 describe generating 
functions for the percolation problem |3^. Similarly, derivatives of the Green functions 



correspond to geometrical correlation functions; for example the spin-spin two point 
function of the Potts model yields the probability that two vertices belong to the same 
cluster. 

We therefore see that the percolation problem is closely associated with a frozen 
model. To obtain the non-trivial model from the frozen model we need to find models 
that are "close to the frozen model." The natural way of doing this is to use an algebraic 
approach. 

2.2.3 Representations of the Temperley-Lieb algebra 

To proceed further it is convenient to think in algebraic terms and turn to a transfer 
matrix formalism |^D[. Consider thus a rectangle of size LxT and propagation in the y 



direction. We define more precisely the geometrical observables in terms of connectivi- 
ties as follows. Number the vertices on a given time slice by 1, 2, . . . , L, running from 
left to right. The space of states 7if^°™ is the set of all possible partitions of 1, 2, . . . , L. 
In a partition of L at time t, two vertices are grouped together if at time t they are 
connected through a cluster that extends into their present and past. For instance, 
for L = 4 the partition (13) (2) (4) correspond to the fact that there is connectivity 
between 1 and 3, while 2 and 4 are not currently connected to any other vertex by 
connections through clusters in their past (see figure 2). Introduce operators, e^, that 
act on For odd subscripts, e2j-i transforms connectivities at time t to connec- 

tivities at time t + 1 by inserting vertical "occupied" bonds at all but the j^^ vertex. 
Thus one has e3|(13)(2)(4) >= |(13)(2)(4) >, and ei|(13)(2)(4) >= |(1)(2)(3)(4) >. 
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1234 
(13)24 



Figure 2: The partition (13)24 means that vertices 1 and 3 are connected 
through a path of occupied edges at previous times (by convention past means 
edges below y = t). Acting with 63 corresponds to adding an occupied vertical 
edge in all but the second column between y = t and y = t + 1 resulting again 
in (13)24. Acting with ci corresponds to adding an occupied vertical edge in 
all but the first column between y = t and y = t + 1 and this destroys the 
connection between vertices 1 and 3, resulting in the partition 1234. 

For even subscripts, e2j, modifies connectivities at the same time, and simply in- 
serts a horizontal edge between the j^^ and (j + 1)**^ vertices. For example, one has 
e2|(13)(2)(4) >=|(123)(4)>. 

These operators satisfy the relations 

e] = \lQej, ejej±iej = ej, [e^, e^] = 0, |j - fc| > 2 (5) 

with ^/Q = 1, and hence furnish a particular representation of the Temperley-Lieb 
algebra |]4ll , 42, 3^. The general form of this algebra can be realized in a number of 



simple models, and, in particular, in the Q-state Potts model. 
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Define Uf"^^^ to be the space of all possible (Potts model) spin states on the L 
vertices across a time slice of the lattice. Introduce the following operators on 7^^°***^: 

(%)<x,<7' = n -^K'^fc) -^(^i'^i+i) (6) 

k 

and ^ 

(e2i-i),,,, = ^ n ^i.^k,(y'k) ■ (7) 

k^j 

The operators Cj satisfy the Temper ly-Lieb algebra for general Q. The transfer matrix 
of the Potts model may be written: 

r = Q'/' n X,, n ^2,-1 (8) 
i=i 0=1 

where 

-I - 1 

^2i-i = + = 1 + ^ . (9) 

The decoupled spin model in the magnetic field is related to this by setting Q = 1 and 
= l+e~^ as in (j^. The transfer matrix of the bond percolation problem is also given 
by (P), but with the Temperley-Lieb representation given at the beginning of this 
subsection. More generally the transfer matrices of the vertex and RSOS models that 
will soon be described, also read as in (|]) and @ but with their own Temperley-Lieb 
representation. 

On a torus, the algebraic structure is slightly more complicated due to the addi- 
tional edge between the first and L^^ vertex. The correct algebra is now the periodic 
Temperley-Lieb algebra. The various possible weights associated to non-contractible 
clusters and the various natural geometrical sectors have precise meaning in terms of 
traces and representation theory . 



2.3 Unfreezing frozen models 

We believe that the foregoing features are generic for an = 2 lattice model. Such a 
model should be closely related to two trivial models: one rigid or frozen, and the other 
some kind of high temperature dual in which the lattice variables are decoupled from 
one another. In between these extremes is the Hilbert space of the N = 2 lattice model. 
Common to all of these models is the representation theory of some underlying lattice 
algebra. The problem of making the N = 2 lattice model is to find the proper observ- 
ables in the decoupled model, or to find the proper "unfreezing" of the frozen model. 
As we will discuss, this may all be thought of in terms of choices of the representations 
of the lattice algebra. We will also find that while the "decoupled" description of the 
model is more intuitive, it is the "unfreezing" process that is easiest to implement more 
generally. To define more precisely the unfreezing procedure, we first want to describe 
two other representations of the Temperley-Lieb algebra, and their frozen limits. 
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2.3.1 Yet more frozen models 



Consider the representation of the Temperley-Lieb algebra provided by restricted solid 
on solid (RSOS) models. (We refer to more general extensions of such models as 
interaction-round-a-face (IRF) models.) These representations have \fQ = 2 cosvr/ (m+ 
1), where m is an integer. The model is defined by introducing heights, £ = 1, . . . , m, on 
the vertices and faces of the original square lattice (figure 3). Represent these heights 
as the nodes of the Dynkin diagram of Am- The transfer matrix now acts on a con- 
figuration space Hfl'^^ whose basis is given by elements = I,. . . ,L} with the 
constraint that neighbouring vertices on the lattice carry heights that are neighbours 
on the Am diagram. Let ve be the components of the Perron-Frobenius eigenvector of 
the Cartan matrix of A^ (or of the incidence matrix of the Dynkin diagram). That is, 
V£ = sin{Ti(!./ {m + 1)). The representation of the Temperly-Lieb algebra is then given 
by: 

The RSOS transfer matrix is then given by using this in (^ and (^. To get the 
representation with (5 = 1, one simply takes m + 1 = 3. The model is rigid, with 
heights alternating between the value 1 on one sublattice, and the value 2 on the other. 
This model thus has two states, depending upon which sublattice takes the value 1 and 
which takes the value 2. 



Figure 3: The RSOS variables are defined on the vertices of the diagonal 
lattice, that is the vertices and the faces of the original square lattice. 

This model has some form of duality with the original non-interaction spin model. 
The non-interacting spin model can be considered as the T —>■ oo limit of an Ising model 



on the medial lattice. The dual [41] is an Ising model defined on the faces and vertices 
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of the original lattice and at T = 0, which is completely "frozen," and has exactly two 
states. 

More importantly for our current purposes, this model has an equivalent vertex 
formulation whose continuum limit is a Gaussian model. This will finally yield more 
direct evidence of the supersymmetry. 

2.3.2 The vertex model and its continuum limits 

We now introduce our final representation of the Temperley-Lieb algebra: the one 
provided by the six- vertex model. Introduce spin variables ± (usually represented by 
arrows) on the edges of the medial graph of the original square lattice (figure 4). Think 
of these spin variables as being basis vectors in C^. The configuration space of the 
transfer matrix is now T-C^l"^^ — {p^^ . To define the operators e^, introduce 4x4 
matrices fi'ij,^;, where i,j,k,l are either + or — , by setting all the entries to zero, except 
the ij,kl entry which is set equal to 1. Then ej is equal to the identity in all but the 
j^^ and (j + 1)*'^ copies of C^, where it is given by 

e, = + - - ^-+,+- . (11) 

The parameter Q of the Potts model is related to the parameter q of the six-vertex 
model by: 

^/Q = q + q~' , (12) 
and so the supersymmetric model corresponds to setting q — exp(i7r/3). 



Figure 4: The six-vertex arrows are defined on the edges of the diagonal 
(medial) lattice. They are also interpreted as domain walls for a SOS model 
whose heights are defined on the vertices and faces of the square lattice (see 
figure 3). 
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The continuum limit of this model is a Gaussian model |45, 19|. More precisely, 



introduce variables on the vertices of the original lattice and in the middle of its faces, 
with values defined recursively by 0„ = 0^ ± 1, that is, the value of above one edge 
of the medial equal to the value under it plus the value of the spin carried by the edge. 
The dynamics of these variables is described by a Gaussian action at large distance, 
and in the conventions where the topological defects are not renormalized, the action is 

A = ^ [ id<pf (13) 



Arc 

with g = 2/3. As a result, the partition function on a torus with doubly periodic 
boundary conditions, 

Z = tr ( r™"^*^" Y , (14) 

is given, in the continuum, by 

Z = Zg(3/2) , (15) 
where Zq is a Gaussian partition function. That is, 

Zg{9) = - E V^-/' , (16) 

e,m(iZ 

where ^ ^ 

hem = 71-^ + rriy/g ) , /^^^ = - | - rriy/g I (17) 



and r]{p) = p^l"^^ n(l ^ p") is the Dedekind eta function, and the elliptic nome is given 
hy p = exp —2nT/ L. Equation (plSj) coincides with the partition function of the first 
minimal (central charge c = 1) = 2 superconformal model with projection on odd 
fermion number. The holomorphic generators of the N = 2 superconformal algebra can 
be written (after normalizing the boson appropriately): 

J{z) = ^90(z) ; Tiz) = - Udm? ■ ^ ' 

If one used the trivial representation of the Temperley-Lieb algebra given by equa- 
tions (I) and d^), instead of the representation given above, it would lead to 

Z = tr(rP°"^)^ = 2^^ . (19) 

In the continuum this partition function is, up to an overall normalization, the partition 
function in the Ramond sector with an insertion of (—1)^, where F is the fermion 
number^. It is therefore constant and corresponds to the Witten index of the model. It 

^In the purely bosonic formulation the fermion number is defined by F = ^'^''(Jo+Jo) ^ 
where J and J are the holomorphic and anti-holomorphic U{1) currents. 
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may also be thought of as the partition function of the continuum topological matter 
model. 

Beside partition functions, various physical observables in the percolation problem 
can be conveniently studied in the N = 2 supersymmetry formalism. Usually the critical 
percolation problem is considered as a c = CFT, i.e. a twisted N = 2 superconformal 
theory. The physical observables of statistical mechanics are then the ones that are 
usually discarded as unphysical from string theory point of view. The probability that 
two edges belong to the same cluster corresponds to an operator with half-integer labels 
in the Kac table of the c = CFT. Only by turning to an = 2 formalism can its 
correlation functions be studied. This operator belongs in the N = 2 formalism to 
the sector "in between" Ramond and Neveu-Schwarz, with spatial boundary conditions 
twisted by (V Its dimension in the twisted theory is h — ^ 



96' 

2.3.3 A practical use of the topological sector 

The identification of a topological sector is not only useful for supersymmetry purposes. 
If supersymmetry is unbroken, the free energy of a non-trivial lattice model, like the 
six- vertex model with q = e^^^^, can be readily computed by turning to the topological 
frozen model. At the special supersymmetric point, it is not necessary to perform a 
Bethe ansatz computation to determine /. This will be discussed further in section 3. 

2.3.4 Quantum group truncation, restriction and freezing 

Because the six-vertex model (with free boundary conditions) has a quantum group 
symmetry, one can use this to reduce the Hilbert space of the model to obtain a new 
truncated model. In the language of height models this correspond to making the RSOS 
restriction. It may also be viewed as the lattice version of the BRST reduction of the 



Gaussian model to the Virasoro minimal models [Q. The first step is to take g to be a 
root of unity. Here we will take q = e*'^/'"+^, for some integer m. One then replaces the 
traces of operators by Markov traces, that is, the trace of an operator O is defined by: 

trM{0) = tr(Og2^^ ) , (20) 

where H is the Cartan subalgebra generator of the Uq{SU{2)) symmetry of the model 



1^, ^ For example, the partition function of the truncated vertex model is defined by 
Z = trM(r^). 

To see how this truncation works, consider the Hilbert space of the vertex model, 
and imagine decomposing it into representations of the Uq{SU{2)) symmetry. The effect 
of the factor of q^^ in (pO) is to weight the contribution of each representation by its 



g-dimension. In particular, this means that only the type 2X representations contribute 
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to the trace^. This means that we can restrict the trace in (20) to a trace over type 
TL representations. The result of this truncation in the continuum hmit is that the 
Gaussian model becomes the minimal model with central charge c = 1 . ^ , . 

° m(m+l) 

To relate the vertex model to a height model one makes a change of basis in the 
vertex model Hilbert space [p!6|, O, H^. At the vertices and faces of the original 



lattice one introduces positive integer heights, in the following manner. One views 
the arrows of the vertex model as defining the basis elements of the spin-i representa- 
tion of Ug{SU{2)). One starts at one side of the lattice with a fixed element of some 
spin- jo representation of Uq{SU{2)); if one tensors this with an element of the spin-| 
representation then one obtains a combination of vectors in the spin- (jo + |) and spin- 
(jo — |) representations. Perform this tensoring successively with the spin-i states on 
the edges across the lattice. After each tensoring with a spin-| state, associate the spin 
of the resulting Uq{SU{2)) representation to the next vertex. At the same time, keep 
track of the total Cartan subalgebra eigenvalue of the state. A sequence of such spins 
across the lattice, along with Cartan subalgebra eigenvalue is a new basis for the vertex 
model Hilbert space (figure 5). The new basis is related to the old by a huge collection 
of Clebsch-Gordon coefficients. The quantum group symmetry means that the transfer 
matrix is independent of the Cartan subalgebra eigenvalue, and so we may discard it. 
The heights, i, are then related to the spins, j , hj i = 2j + 1. 



Figure 5: The geometry used in the vertex-RSOS transformation. Heights 
sit on the vertices of the heavy square. The dotted edges carry the spins of 
the six-vertex model. Heights are viewed as Uq{SU{2)) highest weights and 
neighbours are connected by tensor product with the fundamental represen- 
tation. 

One can rewrite the transfer matrix, and indeed the individual Boltzmann weights, 
in this new basis. They are related to the vertex model Boltzmann weights by quantum 

^The other representations are called type X and are those representations that are 
reducible, but indecomposable. Such representations also have vanishing g-dimension ITt). 
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6j-symbols [^6|, The effect of performing tlie quantum group truncation in tlie 
vertex model is equivalent to restricting the spins, j, of the representations on each 
vertex to those of type II representations. That is, for q = e*'^/™'+^, one restricts the 
heights to the range £ = 1, 2, . . . , m. In making the change of basis one finds that 
because the quantum 6j-symbols vanish at certain strategic values of the spins, the the 
new Boltzmann weights preserve such a truncation. One could naturally ask what would 
happen if one made the change of basis but did not perform the truncation. One would 
find that the some other 6j-symbols would vanish, and as a result some Boltzmann 
weights would become singular. Thus, when g is a root of unity, this particular height 
formulation is necessarily the restricted (RSOS) model. 

For q = e*'^/^, the only type H represent at aions of Uq{SU{2) are spin-0 and spin- 
i. Therefore, as has already been observed, the heights in the RSOS model alternate 
between 1 and 2, and the model is frozen. 

From this perspective it is somewhat clearer what we must do in order to recover 
an = 2 model from a frozen model. We must generalize a process that takes us 
back to the six-vertex model from the RSOS model. Morally speaking, this should 
involve some form of releasing the height restriction, however it is not quite this sim- 
ple since, as we mentionned above, doing this in the RSOS model will result in both 
vanishing and singular Boltzmann weights. The key point is to realize that the solu- 
tion to this problem of singular Boltzmann weights is directly related to the process 
of untwisting the energy-momentum tensor of the continuum topological theory into 
the energy-momentum tensor of the N = 2 superconformal theory. In more concrete 
terms, the fact that the transfer matrix in the six vertex model with free boundary 
conditions commutes with the quantum group means that the transfer matrix and the 



corresponding spin-chain hamiltonian necessarily contain special boundary terms [17 
The continuum limit of this hamiltonian is that of the topological theory. Thus to 
obtain the N = 2 lattice model one has to twist the transfer matrix of the vertex model 



with free boundary conditions so as to remove these boundary terms ||3^. Doing this 
also breaks the quantum group symmetry and thus modifies the spectrum of the spin- 
chain hamiltonian The result is the pure Gaussian model described earlier. An 
alternative way of removing these boundary terms is simply to use periodic boundary 
conditions (this was implicitly used earlier). This works for the six vertex model, but 
not for its generalizations, in which we only want to break part of the quantum group 
symmetry. We will discuss this more extensively in the next section. 

2.4 Where is the supersymmetry on the lattice? 

As stressed above, one of the indications of supersymmetry in any model is the existence 
of a trivial sector. We have also seen that the N = 2 superconformal generators can 
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be explicitly constructed in the continuum limit. One would, however, like to see more 
direct evidence of the supersymmetry on the lattice. For example, one would like to 
be able to exhibit lattice quantities that reproduce the supersymmetry algebra. More 
generally, one would like to find some lattice fermion operators. 



Unfortunately, although some progress has recently been made in the identifi- 
cation of lattice quantities whose continuum limit is the Virasoro algebra, the situation 
is still very unclear for supersymmetry generators. An obvious related difficulty is that 
on the lattice we only see the symmetry Uq{SU{2)) with q = exp(27r/3), while the 
continuum theory is characterized by a pair of quantum groups, the other one being 
Uq'{SU{2)), q' = exp(OT/2), and it is this latter quantum group that has much to do 
with the supersymmetry algebra. For the more general models of next section, we will 
again observe the "wrong" quantum group on the lattice. 

For the c = 1 superconformal model there is, however, an = 2 lattice model based 
on polymers that exhibits a Uqi{SU{2)) symmetry. Although no explicit realization of 
supersymmetry is known in this model, it at least allows a rather satisfying identification 
of lattice fermionic degrees of freedom. Polymers, or self-avoiding walks, are better 
described as the limit of the 0{n) model |^ as n ^ 0. The topological model is simply 
a model with no degrees of freedom at all, and with Z = 1. The non-local geometrical 
observables are obtained by considering properties of self-avoiding mutually avoiding 
walks on the lattice. Typical correlation functions at coupling (3 have the form 

E/?"^^^, (21) 

where the sum is taken over all self-avoiding walks connecting two points (see figure 6), 
and ^Tv, is the number of such walks of length A^. The critical point, where the corre- 
lators decay algebraically, occurs when (3 = Pc, where is the "effective connectivity 
constant," which is defined for large A^ by ~ Pc^^ N >> 1. The proper algebraic 
setting for this model is It has however a Uq{SU{2)) symmetry with q = i 



A simple way of identifying fermionic degrees of freedom is to think of the zero weight 
given to closed loops as being the sum of statistical factors +1, the —1 corresponding to 
the loop carrying a bosonic or fermionic variable. On a torus, the choice of antiperiodic 
boundary conditions for fermions gives a weight 1 + 1 = 2 to some families of loops: 
this allows a simple recovery of the various sectors of the N = 2 theory from the lattice. 
This fermion interpretation allows also a transparent interpretation of the new index of 



14, 15] in polymer terms. 
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Figure 6: The derivative of the spin two-point function in the 0(n) model, 
evaluated at n = 0, becomes a generating function for the self-avoiding walks 
that connect two points. 

2.5 Off-critical N = 2 supersymmetry 

The first, and most obvious property guaranteed by (unbroken) supersymmetry in a field 
theory is that the ground state energy will be exactly zero. This does not immediately 
imply that the bulk free energy of the lattice model is zero since the partition function, 
Z, is ambiguous up to non-universal factors of the form e^'^ , where N is the number 
of "sites" of the model. Equivalently, one can always multiply all of the Boltzmann 
weights by some overall function that is analytic and nowhere vanishing in the regime 
of interest. This will produce the foregoing ambiguity in the partition function. Thus, 
in a supersymmetric lattice model we should expect that the free energy will vanish up 
to logarithms of such analytic functions. 

As we remarked in the introduction, the percolation problem has N — 2 super- 
symmetry in the continuum limit only at its critical point. To have unbroken super- 
symmetry in the scaling region it is necessary that the trivial sector remains trivial in 
the perturbation, and moreover that the free energy per vertex continues to be sector 
independent. If this is so, then the free energy cannot have any singular part. Recall 
that in general the free energy has two parts: 

/ /reg ~l~ /sing ) (^2) 

where /reg is analytic in the variable measuring the distance away from criticality; and 

/sing = A± , (23) 

where A± are amplitudes above and below the critical point, and ^ is the correlation 
length. The singular part of the free energy, /sing, can be identified with Eq/L where 
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Eq is the ground state energy of the quantum theory. If supersymmetry is unbroken 
then one has Eq = 0, and therefore fsmg = 0. 

This lattice criterion can immediately be used to claim that for the polymer problem 
the flow to the dilute region (where /? < /?c in (0)) does not break supersymmetry. On 
the other hand the flow to the dense region (3 > (3c does break supersymmetry. This is 
because for bigger weights, the walks fill a finite fraction of the available space (figure 
6) ^5[, so as long as one loop is allowed there is a non-trivial free energy: = 1 but 
Z ^ e-^'^^ otherwise. Supersymmetry turns out to be broken spontaneously, a fact that 
is possible because of the non-unitarity of the perturbed polymer problem. We will 
discuss this more in a later section. 




Figure 7: If the variable (3 of the generating functions is real and bigger 
than the radius of convergence /3c in infinite volume, the self- avoiding walks 
start to fill the system and properties depend very much on the boundary 
conditions. The free extremities of the walks repel one another algebraically. 



If one starts with the six vertex model at criticality, then the eight vertex model is the 
obvious off-critical generalization. In the regime where the elliptic nome is positive, and 
for q = e*'^/^, the supersymmetry is not spontaneosly broken. (This will be established 
as a consequence of the results in the next section.) Indeed, it has been known for some 
time that the singular part of the free energy vanishes for the eight-vertex model for 
this value of q In the continuum limit, the eight- vertex model can be thought 

of as a theory with Landau- Ginzburg potential W = — XX. In the other regimes, 
where the theory is non-unitary, we expect that supersymmetry will be broken. 

We believe that the foregoing is a generic property of = 2 lattice models: there 
will be an off-critical regime in which supersymmetry is unbroken, and a non-unitary 
regime where its spontaneously broken. 
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2.6 Summary 

At the end of this long, and perhaps somewhat confusing, section we would hke to recall 
the properties that we expect of an = 2 lattice model and summarize them in the 
form of a strategy that can be used in seeking out such lattice models. 

For any Lie algebra, spin and quantum group parameter q there is a natural lattice 
algebra with which solutions of the Yang-Baxter equations are built. This algebra 
has several types of representations, including a RSOS one, a vertex one, and some 
others of mixed type [Q. Choose the value of q such that the RSOS representation is 



trivial. One obtains then a candidate for a topological sector. Study then the complete 
model obtained from the vertex representation or maybe the mixed one. Compute 
in particular its torus partition function and compare it with known N = 2 ones. 
To investigate unbroken supersymmetry away from the critical point, study integrable 
deformations of the critical Boltzmann weights, and (i) check that the free energy has 
no singular part (ii) compute the scaling dimension and U{1) charge of the lattice 
operators corresponding to the integrable perturbation and verify that the continuum 
limit of these operators preserves the supersymmetry. 

3 The coset construction of = 2 lattice models 

Thus far, we have considered only the simplest of the N = 2 lattice models: those whose 
qunatum group structure is Uq{SU{2)) and whose Hilbert space forms a representation 
of the Temperly-Lieb algebra. The critical, continuum limits of these models are the 
N = 2 superconformal minimal models. To go beyond this and obtain theories from 
models with larger quantum group symmetries, and whose continuum limits are the 
N = 2 superconformal coset models, one first finds a formulation of the coset models 
in which the topological sector has a simple lattice formulation. 



Consider the N = 2 superconformal coset model of the form |66 



H 



where G is simply laced andG/H is a hermitian symmetric space. It was observed in 
that there is a natural perturbation of this model leading to an = 2 supersymmetric 
quantum integrable field theory, J^l{G] H). Moreover, the topological subsector of this 
quantum integrable model can be identified with the topological coset conformal field 
theory: 

■ (25) 

In particular, the correlators of the topological sector of Ail{G; H) are all constant, and 
can be written in terms of the structure constants of the fusion algebra of G^- Since 
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lattice analogues of the Gk x Gi/Gk+i are well known, the basic approach should now 
be evident. To see more precisely how the N=2 lattice model is constructed we need to 
elaborate some of the details of the superconformal model. 

3.1 A Coulomb Gas Formulation 

For simplicity, take the level of G in (|^) to be one {i.e. k = Let r be the rank of 
G. Since G and H have the same rank, the representations of the current algebra of Gi 
are finitely decomposable as representations of the current algebra of Hi. Because G/H 
is a symmetric space, the representations of the current algebra of SOi{dim{G / H)) are 
finitely decomposable into representations of the current algebra of Hg^h, where g and 
h are the dual Coxeter numbers of G and H respectively. It follows that Aik{G] H) can 
be thought of as a coset model of the form: 

Hi X Hg_n/Hg_n+i , (26) 

but with a special choice of modular invariant. Because of this equivalence, one can find 
a Coulomb gas formulation that directly generalizes the SU{2) Coulomb gas formulation 
of subsection 2.3.2. For SU{2), this was a simple Gaussian model, here one gets a 
Gaussian model with r free bosons compactified on a scaled version of the weight 
lattice of G. The model has a quantum group symmetry of f/g(i/o) where Hq is the 
semi-simple factor of H. (For hermitian symmetric spaces G/H, the group H has the 
form H = Ho X f/(l), where Ho is semi-simple.) The quantum group parameter, q, is 
the one appropriate to the denominator factor of the coset (^), that is, one takes 

q = e (3-h+i)+h = eiTT . (27) 

This choice fixes the radius of comapctification of the Gaussian model to the "super- 
symmetric radius." 

To reduce the Gaussian model to the requisite coset model one must perform the 
BRST reduction of the field theory, which, on the lattice, is the quantum group trun- 
cation with respect to Ug{Ho). The model knows it origins as M.k{G;H) essentially 
because the Gaussian model is compactified upon the scaled weight lattice of G. The 
field theory also contains two operators, and X~ , that in the Coulomb gas formu- 
lation extend the generators of Ug{Ho) to a twisted form of the affine quantum group 
Uq{G). In the N = 2 superconformal model, these operators are hermitian conjugates 
of each other, they are relevant, and together provide a perturbation that yields a uni- 
tary N = 2 supersymmetric quantum integrable field theory |]68|, 0. This is directly 

^The field theories with higher levels, k, can be obtained by including generalized 
parafermions, and the lattice models with higher values of k can be obtained by fusion. 
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parallel to the situation for non-supersymmetric quantum integrable models obtained 
by conformal perturbation theory: the perturbation is usually one that leads to an 
affine extension of any underlying quantum group structure. Here, however, one needs 
two perturbing operators to make the quantum integrable model, one of which extend 
Uqi^Ho) to Uq{G) and the other extends this to Uq{G). 

The topological twist of any N = 2 supersymmetric field theory can be implemented 
by first replacing the energy momentum tensor by |6^, [70[| : 



= T;:=' + \e,, g^-d^J^ , (28) 

where J^, is the conserved U{1) charge of the N = 2 supersymmetric theory. With 
this energy momentum tensor, two of the supercharges become dimension zero charges 
and can be used as BRST charges to truncate the original Hilbert space down to the 
topological Hilbert space. 

The topological twist shifts the dimension of the operators and X~ so that they 
provide exactly the correct charges to extend Uq^Ho) to Uq{G). The topological energy 
mometum tensor thus commutes with Uq{G). The type H representations of Uq{G) at 
the value of q that we have chosen in (^) are all trivial, and the BRST reduction, or 
quantum group truncation, will result in a rigid topological model. 

In the foregoing discussion we used a minor sleight of hand that we wish to bring 
into the open. In the continuum field theory formulation of a coset model there are 
always two quantum groups, one associated with a numerator factor and the other 
associated with the denominator factor. Thus in ( p6|) there are two quantum groups, 
the one that we have been discussing with q given by (|27|), and one associated with 
the numerator factor of Hg_h with q = e^'^^^ . Two supercharges can be added to the 
latter quantum group extending it to a twisted form of Uq{G), while the extension of 
the denominator quantum group is accomplished by the relevant perturbing operators 
and X~ as described above. To define the topological theory one is supposed to 
use one of the supercharg GS clS 8b BRST charge, and this implicitly suggests that one 
gets the topological theory using the numerator copy of Uq{G). However it should be 
remembered that the BRST reduction of a Coulomb gas formulation of a coset model 
can be done using either one of the two quantum groups, and we have used, and will 
need to use, the denominator copy of the quantum group. 



3.2 Formulating the N = 2 Lattice Models 

To build the N = 2 lattice models one simply reverses the foregoing course. One starts 
either with a vertex model built using the fundamental representation, V, of G, or 
with the IRF model whose heights are the weights of G. The "topologically twisted" 
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transfer matrix is built from the ^-matrix of Uq{G). If one performs the quantum group 
truncation of the vertex model ||T^ [1^ , or the restriction of the height model, using the 
complete Uq{G) with q given by (|27D the result is a rigid lattice model. In the IRF 
formulation each successive height will be the unique level one fusion product of the 
previous height and the representation V. 

To get the N = 2 lattice model one must topologically untwist the foregoing transfer 
matrix, breaking the Uq{G) symmetry to f/g(ifo) in such a way that the two extra 
generators in Uq{G) are given the same spin. One then only performs the quantum 
group truncation with respect to Uq{HQ), or, in the IRF formulation, only restricts 
the heights to be affine highest weights of Hq at level g — h + 1, leaving the U{1) 
direction unrestricted. This is referred to as partial quantum group truncation, or a 



partial RSOS model. This is the procedure as it was originally implemented in |36 



In the Coulomb gas formulation was extensively analysed to precisely define the 
topological untwisting of the transfer matrix, and to provide further evidence that the 
result was indeed an = 2 lattice model. In this review we will use more recent results 
to define the off-critical Boltzmann weights in the IRF formulation of these models 
with G = SU{N). It turns out that in the off-critical formulation the "untwisting" is 
easy to describe and is rather intuitive. We begin by reminding the reader about the 
construction of RSOS models based on the weight lattice of SU{N). 

Introduce an orthonormal basis, fj, j = 1,...,A^, in IR^ and define vectors Cj 
by Gj = /j — -^(/i + . . . + /tv). The vectors ej can be thought of as weights of the 
fundamental of SU (iV). Consider the oriented lattice shown in figure 8. To each vertex 
assign a height of the form: 

Af-l 

A = Ao+Y. ' (29) 

i=i 

where nj G 2 and Aq is an as yet arbitrary vector. This "initial vector," Aq, will play 
a major role in the forthcoming discussion. 

If A is a height at the beginning of an oriented edge and A' is the height at the 
other end of the oriented edge, then we require that A' — A = Cj, for some j. For the 
typical plaquette shown in figure 9, the evolution is defined by the Boltzmann weights 
w{A, A + Ci, A + Cj + Cj, A + ek\u). Introduce the shorthand notation: 

[u] = ^i(7z/|r) (30) 

oo 

= 2plsin{7r-fu) - p")(l - e2"T>")(l - e'^'^^^V) , (31) 

n=l 

where q = e™'^ is the quantum group parameter and p = e^'^*'^ is the elliptic nome. Let 
p denote the Weyl vector of SU{N): 

P = liiN- l)ei + {N- 3)62 + - (iV - l)e^) . (32) 
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Figure 8: A section of the lattice upon which the model is defined. The bold 
zig-zag is the initial time slice, and the arrow indicate orientations of edges. 



The non-vanishing Boltzmann weights are as follows |^ : 



w{A, A + Ci, A + Ci + Cj] A + ei\u) 
w(A, A + Cj, A + Cj + Cj] A + ej\u) 



w{A, A + e^,A + 2ei; A + ei\u) 



[^ + 1] 
[1] 

[{A + p) ■ {a - ej) - 



u] 



[(A + p)-(e,-e,)] 



N / [(A + p) ■ (e. - e,) + 1] [(A + p) ■ (e. - e,) - 1] 
[1] V [(A + p)-(e.-e,P 



) 



(33) 



where i ^ j- 

To obtain the restricted (RSOS) models corresponding to Ge x Gi/Gi+i, one takes 
Ao = 0, 7 = + g + 1) and restricts the heights to the fundamental affine Weyl 
chamber of G^+i- That is, the heights must be affine highest weights of G^+i. If one 
were to make this choice of Aq and 7, but not make this restriction on the heights, 
then the Boltzmann weights would be singular for certain combinations of heights. It 
is elementary to verify that the transfer matrix arising from the Boltzmann weights in 
( |33[ ) preserves this restriction. The topological, or rigid model corresponds to setting 



e = o. 
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Figure 9: In the IRF formulation the heights, A, are associated to vertices 
as shown. 



The lattice analogues ||3^ of the N = 2 superconformal grassmannian models: 

SUi{m + n) X S0i{2mn) 



SUn+iim) X 5f/„+i(n) X f/(l) 



(34) 



are constructed by once again taking the heights as in (^) (with N = n + m). However 
one now performs the partial restriction by requiring that the heights be affine highest 
weights of both of the subgroups SUn+iim) and SUm+i{n) of SU{m + n), but one 
does not restrict the heights in the U{1) direction. The value of q is taken to be the 
"topological value" (|27|) . The topological untwisting is accomplished by choosing 

Ao = — (pg - Ph) t , (35) 
19 



where pc and pn are the Weyl vectors of G and H. For the grassmannian model (|3^) , 
this is simply 

in -\- 71 -\- \ 

^0 = TT/ — r[n{ei + . . . + em) - rn{em+i + ■ ■ ■ em+n)]r . (36) 

2[m + n) 



There are several things to note about this choice. First, the proof in |71| that 
the Boltzmann weights ( P^D satisfy the star-triangle relations works for a general value 
of Aq so this model is indeed still exactly solvable. Secondly, Aq is orthogonal to the 
SU{m) and SU (n) directions and so these sectors of the model behave as if Aq were zero. 
This means that the transfer matrix preserves the partial restriction on the subgroup 
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SU{m) X SU{n). Hence this model is an IRF formulation of a particular form of 
(pGl). Because the U{1) direction is unrestricted, there is, in principle, a danger that 
the Boltzmann weights could be singular for some configuartions. In addition, in the 
physical model, the vector (|35D is purely imaginary, and so one might be concerned 
that the Boltzmann weights may cease to be real in the physical regime. As we will see 
momentarily the Bolzmann weights are still real and non-singular. 

It is convenient to think of the foregoing models from slightly different perspective: 
the heights can be taken on the partially restricted weight lattice of G = SU (m + n) 
with Ao = 0, but the Boltzmann weights must now incorporate the shift (P5|). We will 
henceforth adopt this perspective. The shift by Aq does not modify the Boltzmann 
weights in the Hq = SU{m) x SU{n) direction, but the other Boltzmann weights are, 
significantly modified. The effect of this shift is to replace certain strategic "i^i's by 
and this neatly removes all the potential singularities in the Boltzmann weights owing 
to the unrestricted U{1). In order to describe the results explicitly let [u] be given by 
( |3lD and let {u} is defined by: 

{u} = Ml'^lr) (37) 

00 

= - e'^^^V"^) (1 - e-2"*^V~^) . (38) 

p=i 

Set the parameter 7 to , and take the vector Aq to be zero. The non- vanishing 

Boltzmann weights are then as follows: 

lfl<i,j<morm + l<i,j<m + n then, for i j, one has: 
w{A, A + Ci, A + 2ei; A + ei\u) = """l 



[1] 

[(A + p)-(e,-e,) - u] 
[(A + p)-(e,-e,)] 
w{A,A + ei,A + ei + ej;A + ej\u) = 



w{A, A + d, A + d + Cj] A + ei\u) 



[u] ( [(A + p) ■ (e, - ) + 1] [(A + p) ■ (e, - e,) - 1] 



[1] V [(A + p) ■ (e, - e,)F ; • 

If 1 < i < m or m + 1 < j < m + n, or vice- versa, then one has 

/A A A A I \ {(^ + P) ■ (^i - e,) - u} 

w; A,A + e„A + ei + ej;A + e, M = ./a n i m 

{(A + p) ■ [ei-ej)} 

w{A, A + A + d + Cj] A + ej\u) = 



(39) 



M / {(A + P) ■ (e. - e,) + 1}{(A + p) ■ jd - e,) - 1} ' 
[1] I {(A + p)-(e.-e,)P 



(40) 
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From the Coulomb gas analysis of |^ and as confirmed by the Bethe Ansatz com- 
putations in the foregoing height model at criticality yields the Neveu-Schwarz 
sector of the N = 2 superconformal model. The Ramond sector can be obtained by a 
uniform spectral flow in the U{1) direction. This is easily implemented: one shifts all 
the lattice height appropriately, or equivalently one takes r ^ r + 1 in (|35|) . 



3.3 Some simple properties of these models 

As has already been mentionned, a signal that a model is supersymmetric is that the 
free energy per lattice site is analytic when the supersymmetry is unbroken. For the 
models discussed in the last subsection, this computation is very revealing in that it 
also exhibits and utilizes the connection between the supersymmetric model and its 
topological sector. We will begin, however, by illustrating the foregoing construction 
for G = SU{2), and relate it to the constructions of section 2. 



3.3.1 A simple example 

For G = SU{2) the foregoing lattice model construction collapses to what is basically 
the IRF version of the eight- vertex model [O, The weight vectors ei and 62 of 



the foregoing subsection satisfy ei = —62, and the Weyl vector is given by p = ei. The 
unrestricted heights lie on the weight lattice of SU{2), which we will parametrize by 
an integer, i, where A = Aq + — l)ei. The integer, (£ — 1), is equal to twice the spin 
of the corresponding SU{2) weight. These are exactly the heights discussed in section 
2.3.1. If one takes Aq = and 7 = 1/3 then the quantum group truncated model is 
completely rigid with i alternating between 1 and 2. There are thus two states in this 
model depending upon whether a given site has height 1 or 2. This is the SUi{2)/ SUi{2) 
lattice model. To get the lattice analogue of the N = 2 supersymmetric theory, one 
should only quantum group truncate with respect to the H subgroup, or more precisely, 
with respect to the semi-simple part, Hq, of H. In this instance, this means that one 
performs no quantum group truncation at all, and one therefore has a Gaussian model. 
The height, i, takes the values 0, 1 or 2 modulo 3, and values of 7 and Aq are: 

1 3 

7 = - and Ao = 2 ^ • 

We have thus "unfrozen" the heights i of the RSOS model and have avoided the problem 
of singular Boltzmann weights through the choice of Aq. As described in section 2.3.2, 
this value of 7 means that at criticality and in the continuum the Gaussian model is 
the N = 2 superconformal minimal model, with central charge, c = 1. The off-critical 
N = 2 supersymmetic, quantum integrable model corresponds to the most relevant 
chiral primary perturbation of the conformal model [|T^, and may be thought of as 
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sine-Gordon at the supersymmetric value of the couphng constant. From (^) and (|40|) 
the elhptic Boltzmann weights are the following: 



w{i,i±l,i±2;i±l\u) 
w{i,i±l,i;i±l\u) 

w{e,e±i,e;eTi\u) 



Ml\r) 
Ml\r) 

^i(f|r) / ^4(i(£-l)|r) ^4(1(^+1)1 
Ml\r) [ iMl\r)y 



T 



These Boltzmann weights are precisely the same as those of the A'^^ cyclic solid- 
on-solid models described in In this context the labelling of the model by 
A'^^ refers to the fact that, because of the periodicity of the Boltzmann weights, the 
unrestricted U{1) is, in fact, cyclic and so rather than taking the heights to be in Z one 
can view them as living on the extended Dynkin diagram of A2. 

3.3.2 The free energy 

The easiest way to compute the free energy per unit volume, JF, is to use its analytic and 
inversion properties (see chapter 13 of ^T[). In this review we will only consider regime 
/// of the SU{N) model, i.e. < p = e^^'^ < 1, -\N < Re{u) < 0. (This regime 
corresponds to the unitary, quantum integrable field theory of interest - other regimes 
either have different conformal limits or correspond to perturbations with imaginary 
coupling.) We will also consider the models defined by ( p8| ) - (^) where 7 is now an 
arbitrary, positive, real parameter. 

To compute it is first convenient to multiply the Boltzmann weights (^) - (|40|) 
by a factor of 

2 

g — (n +-) , (42) 

and perform the modular inversion r — After making a gauge transformation 
one then finds that the Boltzmann weights are manifestly periodic under: 

u ^ u + — . 43 

7 

One then writes 

T = - log{K{u)) , 

and requires that log{K) be analytic in a region containing regime III and also be 
periodic under (|43|). Note that one does not impose the other periodicity of the theta 
functions {u m+^) because such shifts of u would go outside the domain of analyticity 
of the free energy. 
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Using the properties of the corner transfer matrix one can show that k{u) must also 
satisfy: 

k{u) k{-u) = h{l - u) h{l + u) , (44) 
k{\ + u) k{X-u) = h{\ - u) h{\ + u) , (45) 
where A = -N/2 for SU{N), and 



Ku) ^ / - ■ (46) 



The function h{u) that appears on the right hand side of these equations depends only 
upon the inversion relations of the elliptic Boltzmann weights. These equations, along 
with analyticity and periodicity, determine log{K{u)) completely. One simply writes 
log{K{u)) as a general Fourier series in e*'^'^"/'^ and uses ( ^4]) and (|45|) to determine the 
coefficients. This is a little tedious but it is straightforward. 

At this point one should note that the foregoing equations and constraints do not 
depend on the choice of Aq. This means that the free energy for all of of the lattice 
analogues of the grassmannian models ( |5^ only depends upon N = m + n, and this 
free energy is exactly that of the topological SUi{N) x SUo{N) / SUi{N) model. This 
will be a general feature of these models, the free energy of the models will be equal 
to that of the rigid topological model, and thus it must be possible to normalize the 
elliptic Boltzmann weights analytically so that the free energy is zero. 

The fact that the free energy is independent of Aq also means that we can use the 
known results for Aq = for the models based on SU{N) |]72[. One therefore has 



where 

sm/i( — (1 — 7)x) sinh{—'~fux) sinh{ '^''^'^~^^ 'jx) 

f{x;u,'y,T) = ^ . X ,.^iN — ^ • (48) 

X sinn[—x) sinh[—^'yx) 



This expression differs slightly from that of ||72[ in that we have subtracted the logarithm 
of the phase (^) so that the result is no longer exactly periodic under (|i5|), but so that 
it does give the free energy for the model defined by the Boltzmann weights in (PPf) and 



To obtain the result as a function of p = e^'^*'^ one needs to perform the modular 
inversion of the second term in (0). This can be done by Poisson resummation (see, 
for example, chapter 10 of HTI], or appendix D of |^). That is, one defines 



r+oc 

/(C) = / e'""'^^ fix) dx (49) 
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and uses the equality: 

fc=oo A;=oo 

E m = E m . (50) 

fc=— oo k=—oo 

The fourier transform of ( ^HD can be performed by closing the contour above or below 
the real axis, depending upon the sign of (, and then summing the residues. This sum 
over residues and the sum in (^) generates an expansion in powers of p. The vanishing 
of sinh{^^'~ix) in denominator of (^T]) gives rise to residues that are proportional to 
piV7 This is the source of the non- analytic behaviour of the free energy as a function of 
p. The residues coming from the vanishing of the other sinh function in the denominator 
of are all proportional to integral powers of p. It is easy to extract the leading 
behaviour as r ^ ioc from this sum over residues. One finds 



4 -2)vr) smC-W) sin{^^ 

sini 



1 



^ — P""" ■ (51) 



1 



For generic values of 7 this means that JF ~ p^^i ^ while from hyperscaling one has 

h. Tl 



JF ~ where ^ is the correlation length. Therefore, we have 



log{0 , (52) 
vr 

as r — * ioo. 

One should note that for 

7"' = Nj + 1, J = 1,2,3,... (53) 

the numerator of (^Sj) vanishes whenever sinh{^^'yx) vanishes. As a result, the Poisson 
resummation of (^7]) gives rise to an expansion in integral powers of p. That is, if 7 
satisfies ( p3|) then the free energy is an analytic function of p. The choice 7"^ = + 1 
corresponds to our supersymmetric models. 

1 

For these special values of 7 it is, of course, no longer true that ~ p^i as r — »• zoo. 
However, continuity in 7 means that the relation (0) is true even at these special values. 

It turns out that for 7 given by (p3D one can easily express the free energy in terms of 
theta functions. Here we will give the result for the supersymmetric model: 7"^ = iV+1. 
The general case is similar. For 7 = 1/(A^ + 1), the second term in ( ^7|) can be written: 

9'Tri °° 1 

— «7(1 - 27) - E ,n ~M - ^~') (^"' - ^"''') ' (54) 
T t^i - p ) 

where p = e"^'^*/'^, z = ^'^'^'^I'^l'^ and w = e^'^^^^'^. This is a standard expansion of 

the logarithm of the ratio of two theta functions^. Finally one can perform modular 

®One simply takes the logarithm of the product formula for the theta functions, expands 
all of the ^017(1 — terms into power series in and then one can perform the 

sum over n to obtain 
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inversions on the theta functions and the final result is 

; . . ^^ ; f Mli^ ' l)\r) \ , f [u - 1] \ 

log{n{u)) = log [-^:^^ j = ^og j . (55) 

Observe that, as promised, this free energy is analytic in the whole of regime ///. 
Indeed, if one multiplies all of the Boltzmann weig hts in (|3|) and §3) by then 
the Boltzmann weights are still analytic in regime ///, and the free energy, JF, vanishes 
identically. 

3.3.3 The off-critical perturbing operators 

For the critical limit of these models one can use Coulomb gas methods to show that 
the continuum limit is = 2 supersymmetric. The free energy computations suggest 
that the off-critical models are also N = 2 supersymmetric in the continuum, but we 
would like to demonstrate this more directly. 

From the analysis of the perturbations of = 2 superconformal coset models, we 
know that the are natural perturbations that lead to A^ = 2 supersymmetric quantum 
integrable models. These operators have the form: 

ijj = G-_,G-,(t) and = G+iG+i^, (56) 

2 2 2 2 

where is a very particular chiral primary field, and is its anti-chiral conjugate. In 
particular, the fields ip and ip have conformal weights and U{1) quantum numbers: 

K - K - - H - 2 + 2(iVTT) 
Qi, = Qi. = -Q^ = -Q^ = 1 - ^^^^ (57) 

We remarked above that in the vertex form of the model at criticality, these operators 
extend Uqi^Ho) to Uq{G). Thus, based on general expectations about lattice models, 
one would expect these operators to be the ones that correspond to the elliptic "de- 
formation" of the critical model. We can see this connection much more explicitly as 
follows. 

Consider the continuous family of models where the initial height vector is taken to 

be 

Ao = —{pG-PH)r, (58) 
19 

where /x is a parameter. When /x = we have the topological model and the lowest 
powers of p in the Boltzmann weights (p3| ) are p° and . If we increase /x, then the 
lowest powers of p change smoothly to p^ and p^^~^^\ Now recall that in the scaling 
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limit, r ioo, we know that r is related to the correlation length by (|52|) , and thus 
p ~ ^ iv+1 . It follows that the elliptic perturbation must involve two operators, and 

_2Njj_ . , 2JV(1-m) 

that their coupling constants must scale as ~ ^ '^+^ and ^ ~ C • The 

corresponding operators therefore have conformal weights: 

When /i = I these two operators have precisely the same conformal dimension and this 
is the conformal dimension given in (|57|). If we now consider the gradual untwisting of 
the energy momentum tensor of the quantum field theory: 

T{z) = T'^'^iz) - fi dJ{z) ; f{z) = f'°P{z) + fi dJ{z) , (60) 

we see that the conformal weights of the operators (^) depend upon the parameter fi 
exactly as in (|59|) . 

In this way one can not only identify the dimensions of the perturbations that lead 
to the elliptic Boltzmann weights, but one can also determine their U{1) charges. More- 
over, the foregoing also demonstrates that smoothly changing the initial lattice vector, 
Aq, causes the dimensions of the perturbing operators to change in exactly the way that 
they should under topological untwisting. As a result one also obtains further confirma- 
tion that one has correctly determined the lattice analogue of topological untwisting. 
A more explicit relation between Aq and perturbed conformal theories can be obtained 



using the method of ||59l . 



4 Scattering matrices 

Since we are considering critical and off-critical N = 2 lattice models, there will of 
course be massless and massive S-matrices associated with the excitations of these 
models. There have been quite a number of papers written on this subject (see, for 
example, |TT|, |TB|, Q), and we will not try perform even a partial survey. 



Our purpose here is to simply make a few remarks about some of the known S-matrices, 
and how their construction fits in with the general strategy of the construction of = 2 
lattice models. 

It is first useful to recall some general observations about conformal coset models and 
their perturbations leading to quantum integrable models. In a generalized Coulomb 
gas description of the conformal model, there are always two quantum groups that 
usually have the same underlying Lie algebra, but have different roots of unity. One 
of these quantum groups is generically associated with the numerator of the coset, 
while the other quantum group is associated with the denominator. The corresponding 
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lattice model only exhibits one of these quantum group symmetries, and it is usually 
the one associated with the denominator of the coset model. The perturbing operators 
that lead to quantum integrable models usually extend this lattice (or denominator) 
quantum group to a larger, affine quantum group. The other (numerator) quantum 
group usually plays a major role in the scattering theory of the quantum integrable 
model. That is, the scattering matrices are usually built out of the i?-matrices of 
this numerator quantum group. The corresponding affine extension of the scattering 
quantum group leads to operators in the theory that commute with the S'-matrix. Such 
non-local conserved charges have been used extensively in the construction of the S- 
matrices. The intuitive reason as to why the scattering theory and lattice quantum 
groups are distinct is that the solitons of the scattering theory must be local with 



respect to the perturbing operators leading to the quantum integrable model fSJ . 

The supersymmetric models also fit this mould. The denominator quantum group 
is indeed that of the lattice model, while numerator quantum group is that of the 
scattering theory. The non-local conserved charges that extend the scattering quantum 
group to the affine quantum group are simply the supersymmetry generators. Thus we 
find that while almost all of our lattice models do not have explicit supersymmetry, the 
corresponding scattering theories do. In other words supersymmetry appears only as a 
dynamical symmetry. 

4.1 Massless S'-matrices and the six- vertex model 



The six- vertex model is diagonalizable by Bethe ansatz and the solutions are clas- 



sified in terms of 1, 2 strings and the anti-string 1 , with a coupling diagram as in figure 
10 (see To study the scattering theory simply observe that for a relativistic QFT 



the thermodynamic Bethe ansatz ||58[ , in addition to the pseudo-particles appearing in 
the solution of the Bethe equations, involves an additional physical particle, as in figure 
10. 

Therefore, if the lattice model has q = exp(z7r/3), the scattering theory has q' = 
exp(i7r/2). From the lattice model, a massless scattering theory can be extracted 
which has the following content (see p9| , |61|]). One has a pair of right moving particles 
which we parametrize as 

eR = PR = —e (61) 
and a pair of left moving particles with 

eL = -PL = ye-^ (62) 
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Figure 10: The first diagram represents the pseudo-particles (solutions of 
the Bethe ansatz) and their couplings at the super symmetric point of the 
lattice vertex model. Recall that for q = e^^^^ one has strings (1,2,. — 1) 
and the anti-string a. In the continuum theory there is an additional node 
for the physical particles, indicated by a heavy dot on the second diagram. 
Therefore the number of Bethe ansatz solutions must differ by one from the 
lattice case, hence the denominator of the quantum group parameter must 
also differ by one from the lattice case. 



where M is a mass scale and ^ is a rapidity. These particles have fermion number equal 
to ±1 and they have factorized scattering with 



'>RR 



m 



( cosh^/2 \ 

isinh^/2 1 

1 isinh^/2 

V cosh^/2 ) 



(63) 



where Z(ff) is a normalization factor ensuring unitarity and crossing symmetry, while 
the LR scattering is trivial 

Slr = Srj^ = \ (64) 
Acting on the right particles h, d we have 



Q-\u{d) >= ^M/2e^l^\d{e) > 
Q+\d{e) >= JM/2e^/^\u{9) > 



(65) 
(66) 

and same type of relations with left particles and generators. The supersymmetry 
algebra is 



Ql = ^Ql^ e = {Q+, Q-} = {Q-, Q+} = 

{Q+,Q+} = {g_,g_} = o 

{g+,g_} = Me^ Q_} = 2Me-^ 



(67) 
(68) 
(69) 
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4.2 Massive 5- matrices 



There are several ways in which one can extract the S'-matrices of the massive quantum 
integrable models that appear in the continuum limits of the N = 2 lattice models. We 



will briefly discuss two of them here. The first approach is that taken in |Tl], [T^ , and 
the starting point is to use the Landau- Ginzburg description of the model to determine 
the ground-state and soliton structure. Since the solitons are generically doublets of 
the superalgebra, one can then use the Bogomolnyi bounds in the superalgebra to 
determine the soliton masses. The solitons also have fractional fermion numbers that 
can be determined from the Landau-Ginzburg potential. This information, combined 
with the fact that the ^-matrix must commute with the supersymmetry, be unitary, 
satisfy the Yang-Baxter equations, and obey crossing symmetry is more than enough 
to determine the S'-matrix completely. This has been done for a large class of models 



12], but there is an even larger class for N = 2 quantum integrable models for which 



the S-matrices are not known. The first reason for this is that one can easily show 
that in this larger class the solitons alone cannot form a kinematically closed scattering 
theory [^, This almost certainly means that breather states must be included 



in the spectrum. This does not pose, a priori, any obstacle to the construction of 
the S'-matrices, but as we will describe below, there could be some further conceptual 
problems to be solved if one is to determine the S'-matrices for these more general 
models. 

We will describe the second method of obtaining the S-matrices in more detail here 
since it is rather more in the spirit of our approach to = 2 lattice models. The 
basic idea is to once again use topological models and a modification of quantum group 
truncation. 

The scattering matrices are well known for the quantum integrable models obtained 
from the natural perturbations of the Gk x G^/Gk^i coset models (see, for example, 
[0, ^ ^ ^ |79|, ^). For simplicity we will only consider the models with with 
G = SU{N) and the level k equal to one. The S'-matrices can then be obtained from 
RSOS restrictions of affine Toda S'-matrices. Let Wj, j = 1, . . . , N — 1, denote the 
fundamental representations of SU{N) {i.e. the anti-symmetric tensors of rank j). 
Each of these representations defines a multiplet of solitons of affine SU (A^)-Toda, and 
these solitons can be labelled by the weights of the representations. 

Define Sjf\6,q) to be the two-body S'-matrix for the scattering of particles in Wi 
with those in Wj: 

S^f\e, q): W,® W, ^ W, ® Wi , (70) 

where, as usual, 6 is the relative rapidity of the solitons and q is the quantum group 
parameter. Requiring the S'-matrix to commute with the Uq{SU{N)) symmetry leads 
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to the following result: 



Sf\e,q) = Me,q) R,,{e,q) . (71) 



The last term, Rji{9, q), is the standard i?-matrix for the quantum group Ug{SU{N)) in 
the principal gradation^. The spectral parameter, x, and the quantum group parameter, 
q, of the i?-matrix are related to the rapidity, 9 and the Toda couping constant, /?, by: 

f fA , 

q = exp I I ) ^ = exp I - I where 7 = ^ _ . (72) 



The scalar factor, Vji{6,q), in ([7T|) is the minimal factor that makes the product vjiRji 
crossing symmetric and unitary. This is relatively easy to compute and can be found 
in |]7S|, |7^. The additional scalar factor Xji{6) is a CDD factor, and it satisfies crossing 
and unitarity by itself. This factor contains all of the necessary poles for closure of the 
bootstrap with the spectrum of soliton masses. It turns out that this CDD factor does 
not depend upon the value of the Toda coupling and so can be evaluated by going to 
the limit {jS 1) in which the Toda model can be related to a Gross-Neveu model |7S 



To get the model corresponding to the perturbed coset model one must first incor- 
porate the background charge, which is done in an analogous manner to the twist in 
the Markov trace of subsection 2.3.4. Let H be some subgroup of G = SU{N), where 
H also has rank — 1, and let pn be the Weyl vector of H. Define Sjf^^\6, q) by 

Sf'"\OA) = [x~P"-'' ^ x~^"-'') Sf^ (^x^"-^' ® xP"-^') , (73) 

where h represents the vector of Cartan subalgebra generators of the quantum group. 
This conjugation only affects the i?-matrix part of the ^-matrix, and for H = G, it 
converts the principal gradation to the homogeneous gradation. 

By construction, the S-matrices Sjf^^\6,q) commute with the action of the finite 
quantum group Uq{H). These generators act in the more familiar rapidity independent 
fashion. Therefore one can use the Uq{H) symmetry to restrict the model exactly as in 
the lattice model. 

For each weight, fx, of Wj we introduce a formal operator Kj^\9). These operators 
satisfy an S'-matrix exchange relation: 



^The i?-matrices that one usually encounters in the literature are written in the homo- 
geneous gradation, which is related to the principal gradation by a trivial automorphism, 
which will soon be described. 
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Let T denote the multi-particle fock space generated by the formal action of the oper- 
ators Kjf\9) on the vacuum. The space is an Uq{H) module, and reducible (for q 
not a root of unity) : 

^ = e V^'^'^ , (75) 



where V^^°- ^ is an Uq{H) module of highest weight A^^-*. Since the Kl3\e) act on J^, 
one can consider their reduction 



K%..m-- V^'^'^^V^'^). (76) 
These operators satisfy the exchange relation: 

K%,„M A ,„,{e,) = (77) 



AH) 



E [sf"'\d.<l))LL K%),Ue,) K%^,„,{9,) . (78) 



The S-matrix for the kinks in this equation is the SOS form, and the foregoing con- 
struction is once again the vertex/SOS correspondence. 

As before, the restriction amounts to taking g to be a root of unity and imposing 
a limitation on the allowed highest weight labels A*-"^-*. To get the S-matrix for the 
preturbed coset model Gi x Hi/Hi^i one must take the Toda coupling so that = 
{£ + h)/{i + h + 1), where h is the dual Coxeter number of H. Therefore, one has: 

q = - expi-in/ii + h)) . (79) 

Note that this is the root of unity appropriate to the quantum group structure of the 
numerator quantum group of the coset modeQ 

To get the N = 2 supersymmetric models one chooses H so that G/H is hermitian 
symmetric, and takes i = g — h. Thus one has q = — e~*'^/^. If one were to take 
H = G one would once again get the topological coset model, but with the foregoing 
choice of q, the only Uq{G) representation with non- vanishing g-dimension is the singlet 
representation. There are thus no solitons in the topological model. Once again one 
finds the N = 2 model by modifying the quantum group truncation procedure. One 
can also verify that under the twisting process (|73D , the two generators that extend 

^"^In the earlier sections of this review, where we were discussing lattice models, we took 
q — e'^/^^, where M was some positive integer. We could equally well have taken q to 
be defined hy q = — e^™/*^. Here, however, the affine extension of the quantum group is 
now playing a more central role and this means that one has to make a specfic choice for 
q. This choice is, of course, convention dependent and we are employing the conventions 



of y,|75| 
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Uq{H) to Uq{G) have a rapidity dependence that makes them spin-i charges. If the 
quantum group truncation procedure parallels that of the conformal theory, then these 
operators will become the supersymmetry generators of the restricted model. 

This process works beautifully and simply for G = SU{N) and H = SU{N — 1) x 
f/(l). For this choice, the fundamental reprexentations Wj decompose into a direct 
sum of two representations of Uq{SU{N — 1)). In the RSOS truncation using Uq{H) 
one therefore finds that each Wj yields two solitons, and these form a supersymmet- 
ric doublet. The U{1) charge becomes the (fractional) fermion number, and the two 
generators that extend Uq{H) to Uq{G) act as the supersymmetry on this doublet. We 
have thus incorporated the supersymmetry directly into the affine quantum group. The 
^'-matrices for the N = 2 model can then be easily extracted from the Toda S'-matrices 

foi G = su{N) ig. 

The surprise comes once one tries this procedure for more general models of the form 
(H), for example, for G = SU{A), H = SU{2) x SU{2) x f/(l). The problem is that the 
quantum group truncation does not respect canonical supermultiplet structure. In these 
more general models all of the Wj^s give rise to RSOS solitons with fractional fermion 
numbers that agree with the Landau- Ginzburg picture. All of the RSOS solitons coming 
from one of the W/s are, of course, transformed into each other by the generators 
of Uq{G). The two generators that extend Uq{H) to Uq{G) are indeed spin-| and 
commute with the S'-matrix. The problem is that there does not seem to be an RSOS 
truncation upon which the algebra of these generators becomes that of the standard 
N = 2 superalgebra. For example, some Wj^s would give rise to a three dimensional 
"supermultiplet." (This happens in the SU{4) example above for the six-dimensional 
representation of S'f/(4).) 

There are several possible resolutions of this issue. It is quite conceivable that the 
Toda approach breaks down, or it might be that there is some supersymmetry anomaly. 
It is also possible that the integrable model gives rise to two ^-matrices, one in which 
the supersymmetry acts locally, and one in which it does not. It is certainly of interest 
to resolve this problem. As first sight, the possibility of a supersymmetry anomaly 
seems the least likely explanation. However, such an explanation is 

perhaps not quite as outrageous as it sounds, after all, the supersymmetry would 
still be a symmetry of the S'-matrix, but it would simply have non-trivial braiding 
relations with single soliton states. Such a phenomenon is already present in the models 
that we do understand: the fractional fermion number of single solitons means that the 
supersymmetry action on single solitons must involve extra phases. It is just conceivable 
that the more exotic models exhibit a more involved (perhaps non-abelian) form of this. 
At any rate, it is highly desirable to compute the S'-matrices for these models. It is also, 
perhaps, significant that these "exotic" models are also precisely the ones for which the 
S*- matrices are not yet known because the solitons cannot, by themselves, form a closed 
scattering theory. 
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5 Other issues 



5.1 Spontaneous breaking of = 2 supersymmetry 

In lattice models there is a natural direction of perturbation that does not break su- 
persymmetry explicitly, but does so spontaneously. This is made possible by the non- 
unitarity of the perturbation. For example, consider the polymer problem discussed 
earlier. The critical Landau-Ginzburg potential in this theory is W = X^, and the 
perturbation that changes the weight (3 of monomers in the generating functions corre- 
sponds exactly to adding a multiple of X to this potential. However if one writes out 
the F-term term in the supersymmetric action, and properly incorporates the coupling 
constant, one obtains: 

Sf = const. (^J £z(fe (x^^ + - f^Y'^ x) + (80) 



j £z(fe (X^ + - f3f/^ X) . (81) 



If /3 < /3c one gets the usual physics of the off-critical quantum integrable model with 
superpotential W = X^ + X. On the other hand, if P > Pc one gets very different 
physics because the perturbation is purely imaginary, and the model is non-unitary. 
This is easily analyzed qualitatively. For P < Pc polymers are "very tiny" (their fractal 
dimension is less than two) so the free energy calculated in the Ramond sector, where no 
polymer at all is allowed, is the same as the free energy calculated in any other sector, 
where some non-contractible polymers are allowed. For P > P^ the partition function 
of the Ramond sector is still Zji = 1 but as soon as a polymer is allowed it fills most of 
the available space and Z grows exponentially with the size of the system. This means 
that there are infinitely many level crossings between the ultra-violet and infra-red fixed 
points in the Ramond sector, and although the states of vanishing energy always stay 
there, they are infinitely far from the true ground state of negative energy in the infra- 
red. In the infra-red the bosonic degrees of freedom of the theory have become massive 
and the fermionic ones are a simple (topologically twisted) Dirac fermion, describing 
the physics of "dense polymers" . 

Interestingly, the position of the level crossings can be related to the poles of the 
Painleve III differential equation Consider the system on a cylinder of radius R 



and length T. The renormalization group variable z is proportional to {Pc — P)R^^^. 
Recall that the generalized supersymmetry index is defined by |T5| : 

Q = 4tr (e-^^ F{-r ) (82) 
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where, from the point of view of the hamiltonian evolution, the "time" is in the R 
direction now. This index can be written: 

where u is the solution of the Painleve III equation: 

(fu 1 du 9 ^ 



with specific asymptotic behaviour |15|. One can show that the level crossings occur at 



the poles of Q on the negative z axis. One can also show that these poles are simple. 



and that there are an infinite number of them |]63 |. 



5.2 Ultra-high temperature Hmits: the non-interacting models 
revisited 

The trivial non-interacting spin model of section 2 has even more to it than was made 
evident there. Recall that the spins sit on the edges of the original square lattice (figure 
1). Suppose that we now draw, on the diagonal lattice of figure 3, 
occupied edges between spins of opposite value. 

These edges correspond to a Bloch wall, and the loops created by them can traverse 
each edge of the lattice at most once. (See figure 11.) These are the usual contours of the 
low temperature Ising model, but considered in the high-temperature phase (actually, 
in the limit of infinite temperature). We now consider the geometrical properties of 
these loops. First it should be clear from preceding discussion that they are identical to 
the properties of the boundaries of percolation clusters (or hulls). The most interesting 
quantity is the fractal dimension of these loops, which is related to the algebraic decay 
of the probability that two edges belong to the same loop. That is, if the latter decays 
as r~^^ then the fractal dimension of the loop is Df = 2 — x. One can show that the 
exponent is given by x = 2h, where h the conformal weight, in the twisted theory, of 
the Ramond ground state with the lowest charge. That is: 

H = i-if-i) = i (85) 

24 2 V 6/ 8 ^ ^ 

and thus Dj = j. This is midway between the fractal dimension for brownian motion 
(Df = 2) and the fractal dimension of self-avoiding random walks (D/ = |). It is 
interesting to observe that, from the point of view of the low-temperature expansion 
for the Ising model, another geometrical quantity is very natural: the probability that 
two edges are extremities of the same open line. This corresponds to the two-point 
correlation function of disorder operators. This function goes to a constant at large 
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Figure 11: The Bloch walls for Ising model at high temperature. Filled dots 
represent states of spin +1 and empty dots represent states of spin —1. The 
walls separate domains of opposite spin. 

distance, and so the conformal weight is zero. This high-temperature limit is therefore 
rather interesting because it rather naturally leads to the existence of a non-trivial 
operator of vanishing dimension in a "lattice topological sector." 

Now we have a natural direction to study: the properties of Bloch walls in high 
temperature phases of spin models. The topological sector corresponds to the original 
spin variables. The non-topological one to the geometry of the Bloch walls. As another 
example consider the three state Potts model. There are three states per vertex, repre- 
sented by empty dot, full dot or square in figure 12. The Bloch walls carry one arrow 
if, when one faces along the arrow, there is either (i) an empty dot on the right and a 
full dot on the left, or (ii) there is a square on the right and an empty dot on the left. 
The Bloch walls carry two arrows if there is a square is on the right and a full dot on 
the left of the arrows. With these rules there is obviously current conservation. 

By analogy with the foregoing Ising model, there are two natural geometrical cor- 
relation functions that are constant at large distance: the probabilities that two edges 
are extremities of the same open line carrying either one or two arrows. Similarly there 
are two natural correlation functions that decay algebraically: the probabilities that 
two edges belong to the same loop carrying one or two arrows. In this example there 
is charge conjugation symmetry that in fact relates simple and double arrows so these 
pairs of correlation functions are in fact identical. The obvious topological model be- 
hind this is the SUi{3) / SUi{3) model. The fractal dimension of the contours can be 
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Figure 12: The Bloch walls for the three-state Potts model at high tem- 
perature. By convention the full dot, empty dot and square correspond to 
a — 1,2,3 respectively. The walls carry one or two arrows depending on the 
difference of spins between both sides. By convention the highest spin is to 
the right of the arrow(s). 

reasonably conjectured to be given by the formula similar to that for the Ising model, 
leading to = y. 

More generally one can conjecture similar properties for the contours in the high 
temperature phase of the models and SUi{n + l)/SUi{n + 1) theories. The con- 
formal weight h will he h — ^~^(^|) = §> where c = is the central charge of the 
-n}^ N — 2 superconformal minimal model. This leads to a predicted fractal dimension 

6 Conclusions and open problems 

The subject oi N — 2 lattice models has, as yet, not produced the same kind of nice 
structures as in the continuum N = 2 field theories. In particular it is important 
to emphasize that there is no real algebraic simplification for the computation of most 
physical properties. This suggests that, despite the large number of models exhibited so 
far, some better ones still remain to be discovered. In addition, very little is understood 
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about the "lattice supersymmetry" and the way these N = 2 lattice models provide 
discretizations of field theories with fermions. A step in that direction would be to 
identify supercharges on the lattice and study their anti-commutator in the scaling 
limit. This might be accomplished in the way that has recently been used for the 
Virasoro algebra [^6] . 

On the other hand the study of = 2 lattice models has already provided many 
interesting physical ideas, in particular concerning the topological models. We believe 
that N = 2 theories are the natural framework to describe most geometrical statistical 
mechanics models and more generally the properties of Bloch walls in non-critical spin 
models. It is also amusing to notice that what is "physical" in field theory and condensed 
matter is rather different. For instance spontaneous breaking of = 2 supersymmetry 
seems to be a very natural feature in lattice models. It is rather remarkable that the 
associated level crossings are encoded in the poles of differential equations of Painleve 
type. 

While finding a lattice version of supersymmetry is probably the most important 
issue to resolve, there are many other open problems in the subject of A^ = 2 lattice 
models. One can try to establish our observations about Bloch walls in generality, or find 
evidence to support these observations through numerical simulations. As mentionned 
in section 4, there are unresolved issues about supersymmetry and the scatering matrices 
in the more complex N = 2 quantum integrable models. (Some suggestions about 
these scattering matrices have recently been made 0].) Simple corner transfer matrix 
computations have been done for A^ = 2 lattice models |^8|, ^ . The results yield N = 2 
superconformal characters as one would expect, but the computation was done by the 



traditional method |^ . It would be interesting to see if the remarkable properties of 
N = 2 supersymmetry could be exploited directly to reach this result in a much simpler 
manner. 

There is also the question of more general lattice correlation functions. In the 
continuum limit, at criticality, one can obtain Knizhnik-Zamolodchikov equations for 
the correlation functions in A^ = 2 superconformal coset models. At first sight, these 
equations are no simpler than those for non-supersymmetric coset models. However, 
once one puts in the special parameter values for the correlators of chiral primary 
fields one finds that one has saturated some bound and the solution of the differential 
equation is particularly simple. For example, if the equation is hypergeometric, the 
correlator of chiral primaries is usually one of the special rational solutions. It would 
be most interesting to see if similar simplifications occur in the g-deformed Knizhnik- 
Zamolodchikov equations for the lattice correlators of p3[ . 

Much of the motivation for these open problems is simply to find some N = 2 
supersymmetry miracles in exactly solvable lattice models. We therefore leave it as 
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a general challenge to find a another formulation of these models, perhaps using the 
Landau- Ginzburg structure more directly, in which the N — 2 supersymmetry leads to 
dramatic simplifications. 
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